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Abstract: The phenomenology and prospects for a discovery of R-parity violating 
Supersymmetry at HERA is analysed. Emphasis is put on the direct resonant 
production of squarks by electron-quark fusion and all possible subsequent decay 



> 

■^j- . modes of the squarks are considered. In particular, the full consequences of the 

mixing in the supersymmetric gaugino-higgsino sector are taken into account. A 

q \ rich phenomenology emerges for HERA which offers a unique sensitivity to new 

R-parity violating couplings and good discriminating power against free parameters 
of the theory. 

1 Introduction 

Supersymmetry (SUSY) which fundamentally links fermions and bosons is likely to be chosen 
as an essential property of a true theory beyond the Standard Model (SM). Among the most 
compelling arguments for a SUSY world are the fact that local supersymmetric transformations 
are fundamentally related to generators of space-time translation (hence necessarily incorpo- 
rates gravity) [Q , the possibility to "explain" the hierarchy between the electroweak mass scale 
and the Grand Unification or Planck mass scale, and the stability of a softly broken SUSY 
which "naturally" avoids the arbitrary fine tuning of the parameters which is necessary in the 
SMi 

A natural framework for SUSY searches is provided by the Minimal Supersymmetric ex- 
tension of the Standard Model (MSSM) |J which has predictive power within a finite and 
well defined set of free parameters and has neither been proven nor falsified by experimental 
observations. The latter is a non-trivial status given the remarkable precision tests of the SM 
at the LEP collider over the recent years. It might have to do with the fact that quantum cor- 
rections due to the sparticles which otherwise respect all gauge symmetries of the SM tend to 
be small rendering indirect observations difficult. It is in addition possible that direct searches 
for particles of the minimal field representation offered by the MSSM have at least partly failed 
because they were looking at the wrong phenomenology. 

The most general Yukawa couplings in a supersymmetric theory which is gauge invariant 
and minimal in terms of field content can be written M in the compact formalism of the 
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superpotential as Wsusy = Wmssm + H 7 ^ • The Wmssm contains terms which are responsible 
for the Yukawa couplings of the Higgs fields to ordinary fermions. The W^ is given by: 

w % = /-,/•-•/.. + K jk LiQjD k + ,\y ; /y; A . (i) 

where ijk are generation indices of the superfields L, Q, E, D and U. The L and Q are left- 
handed doublets while E, D and U are right-handed singlet superfields for charged leptons, 
down and up-type quarks, respectively. The A and A' terms induce lepton number violation 
while the A" terms induce baryon number violation. In the strict MSSM framework, one imposes 
that the SUSY theory be also minimal in terms of allowed couplings and all terms are 

avoided by imposing a strict conservation of the R-parity defined as R p = (—i s ) 3B + L + 2S = \ 
(for particles) = —1 (for sparticles), where S denotes the spin, B the baryon number and L the 
lepton number. Imposing this discrete symmetry is a somewhat ad hoc prescription. Another 
viable || and less restrictive discrete symmetry is the 5-parity which imposes only baryon 
number conservation (i.e. A" = 0). Vanishing A" couplings is sufficient to avoid unacceptable 
n — n oscillations and fast proton decay. Moreover from the cosmological point of view, the 
observed matter/antimatter asymmetry imposes much more severe constraints on A" than on A 
or A' ||. Finally, 5-parity is favoured over i?-parity conservation in a large class of superstring 
inspired models ||. It is also interesting to note that the A and A' terms in ([[J) which have 
no equivalent in the SM arise in a fundamental way from the fact that S'f/(2)-doublet lepton 
superfields have the same gauge quantum numbers as the Higgs supermultiplets. 

The ep collider HERA which provides both leptonic and baryonic quantum numbers in 
the initial state is ideally suited for $ p searches. This was realized long ago and was first 
investigated theoretically in the context of the previous HERA Workshop which motivated 
early experimental searches |J. The cases A' 7^ which could lead to resonant production of 
squarks via e-q fusion offers of course the most exciting prospects. Recent investigations || [K], 
|TT| 1 have shown that a new and rich phenomenology (different for e~ and e + beams) emerges 
when considering the full complexity of the mixing in the gaugino-higgsino sector of the theory. 
This is studied in more details in this contribution in view of future high luminosity runs at 
HERA. 

The case of associated e-q production at HERA followed by the $ p -decay of the sfermions 
has already been studied in detail and also in view of high luminosity runs at HERA in W2 



Via this process one can probe significantly smaller Yukawa couplings than via the resonant 
production but only at smaller sfermion masses. 



2 Phenomenology of ftp SUSY 
2.1 Modelling and Free Parameters 

The X'ijfrLiQjDk terms in the $ p extension of the MSSM correspond in expanded field notation 
to the Lagrangian 



C'LiQjDk — Kjk 



+W k R + u L d{d k R + {vlTVn +h.c. (2) 
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where the superscripts c denote the charge conjugate spinors and the * the complex conjugate 
of scalar fields. Among the 27 possible A - jfc couplings, the cases i = 1 can lead to direct squark 
resonant production and are thus of special interest at HERA. These cases are studied first in 
this paper assuming conservatively that one of the A' dominates. 

The masses of the scalar quarks and scalar leptons, bosonic sparticle partners of the SM 
fermions, are treated here as free parameters. In the gaugino-higgsino sector, there are four 
neutralinos Xi = 1 • • • 4) which are mixed states of the photino 7, the zino Z and the 
super symmetric partners H® and H® of the two neutral Higgs fields. Two charginos xf (j — 1> 2) 
are mixed states of the winos and of the SUSY partners of the charged Higgs fields. The 
masses and couplings of the x° an d X ± are calculated in terms of the MSSM basic parameters : 

• Mi and M2, the U(l) and SU(2) soft-breaking gaugino mass terms; 

• /i, the Higgs mixing parameter; 

• tan /3, the ratio of the vacuum expectation values of the two neutral Higgs fields. 

The number of free parameters is reduced by assuming a relation at the Grand Unification 
(GUT) scale between Mi and M 2 (see Appendix for detail). No other GUT relations are used 
and in particular the gluino (g) mass is left free. 

We moreover consider the following simplifying assumptions : 

• all squarks (except the stop) are quasi-degenerate in mass; 

• the lightest supersymmetric particle (LSP) is the lightest neutralino x\'-> 

• gluinos are heavier than the squarks such that decays q — > q + g are kinematically forbid- 
den. 

It should be made clear that there are no compelling cosmological constraints in $ p models 
which impose that the LSP be neutral and colourless. Other possible choices for the LSP 
(e.g. g or x^) would not significantly change the search and analysis strategy and will only be 
briefly discussed. In l/t p models, in contrast to the strict MSSM, the LSP is generally unstable. 
This leads to event topologies which differ strongly from the characteristic "missing energy" 
signal due to LSP's escaping detection in the MSSM. Hence, except for exclusion limits derived 
from indirect searches (e.g. from the intrinsic width of the Z°), the mass constraints obtained 
in the MSSM framework do not apply directly in ty p models. The search for T^ v squarks is 
"complementary" (hence mandatory) to that performed in the strict MSSM framework. 

2.2 Squark Production 

The resonant squark production mode through direct e-q fusion is illustrated in Fig. [I] for 
A' m 7^ 0. By gauge symmetry, only -u^-like or rf/j-like squarks (or their charge conjugates) can 
be produced in ep collisions. The production of "left" squarks (i.e. supersymmetric partners 
of left-handed quarks) is the dominating process if HERA delivers positrons, since the fusion 
occurs via a d valence quark. On the contrary, with electrons in the initial state, mainly 
"right" squarks are produced. This dichotomy has important consequences since "left" and 
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e + (e + , v e ) e + u,d (d) 



A'n\ — _h_/'X 111 




d (u)/ (d R ) \ d (u, d) d (u)/ (d R ) \ x °, x + ( x °) 
(a) (b) 



Figure 1: resonant production of Ul or d R squarks in e + p collisions with subsequent (a) $ p 
decay or (b) gauge decay involving a (generally) unstable gaugino-higgsino (x° orx + )- 

"right" squarks have different allowed or dominant decay modes as will be seen in the following 
sections. In particular, new exotic final state topologies might have sizeable contributions in 
e + p collisions. 

HERA offers a high sensitivity to any of the nine A' ljfc couplings, in contrast to most indirect 
processes. The production processes allowed for each A' lj7c are listed in Table p] for an e + beam. 
Squark production via X' l j l is especially interesting in e + p collisions as it involves a valence d 
quark, whilst X' nk are best probed with an e~ beam since squark production then involves a 
valence u quark. 

Figure H shows the production cross-sections via A' m for ul and d^, and for cl via A' 121 , each 
plotted for coupling values of A' = 0.1. In the narrow width approximation, these cross-sections 
are simply expressed as 

7T w9 ,,M 2 , 



where y/s = ^jAE^E^ ~ 300 GeV is the energy available in the CM frame for incident beam 
energies of E® = 27.5 GeV and E® = 820 GeV, and q'{x) is the probability to find the relevant 
quark (e.g. the d for Ul or cl and the u for d^} with momentum fraction x = M 2 / s ~ M~/s in 



Table 1: Squark production processes at 
HERA (e + beam) via a R-parity violating 
A ljfc coupling. 
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the proton. Hence the production cross-section approximately scales in A' 2 . The full kinematic 
domain can be probed at HERA for couplings weaker than the electromagnetic coupling (i.e. 
A 2 /47r < a em ) given an integrated luminosity of ~ 500pb _1 . 
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Figure 2: Squark production cross- 
sections in e + p collisions for a coupling 
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2.3 Squark Decays 

The squarks decay either via their Yukawa coupling into ordinary matter fermions, or in a 
first step via their gauge coupling into a quark and a neutralino xl (i = 1 ... 4) or a chargino 
Xj (j = 1,2). The former modes are henceforward called "squark ft p decays" and the latter 
"squark gauge decays" . 

ftp decays of squarks: 

In cases where both production and decay occurs through a A' ljfc coupling (e.g. Fig. |TJa for 



A' m 7^ 0), the squarks behave as scalar leptoquarks Jl3], [L4|. For A' m ^ 0, the oIr resembles on 
event-by-event the S° leptoquark and decays in either e + + u or v e + d while the ul resembles 
the 5*1/2 and only decays into e + d. The partial decay width reads : 

T q^ft p = T u L ^e+d = r 4^ e+s = T g R -, u d = Tq^ X '? 1iM Q ^ 

so that squark ft p decays will mainly contribute at high mass for large Yukawa coupling values 
A'. Hence, the final state signatures consist of a lepton and a jet and are, event-by-event, 
indistinguishable from the SM neutral (NC) and charged current (CC) deep inelastic scattering 
(DIS). 

Gauge decays of squarks: 

The MSSM Lagrangian contains terms coupling a sfermion to an ordinary fermion and a 
gaugino-higgsino. The partial widths for squark gauge decays depend on MSSM parameters 
via the composition of the neutralinos or charginos. 
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Both q L and qn squarks can decay via q — > The partial width of the q — > qx\ decay is 
calculated to be 

1 (M 2 " 

1 ' l2 , D 2\i# 1 Xi 



where : 



A = 9 Jf qN "» > i? = eeXi+^(0.5-e g sin 2 M^, (6) 

and where (-/V?-) is the transport matrix which diagonalizes the neutralino mass matrix (see 
Appendix for detail) in the A — W 3 (7 — Z basis). In practice, the Xi masses and the exact 
values of the "chiral" couplings A and B depend on the relative 7, Z and H components of the 

xl 

The dependence of the Xi mass on the \x parameter is shown in Fig. Qa for fixed M2 and 
tan p. The dominant component (7, Z or H) of the lightest state Xi is shown in Fig. Hb. More 
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Figure 3: (a) Physical masses of the x° Q-nd xt a function of // for tan/3 = 1 and M2 = 
60 GeV; (b) Main component of the LSP for tan/? = 1 



details on the way the nature (and masses) of the various neutralino states depend on the basic 
MSSM parameters M 2 , [i and tan/3 can be found in the Appendix. 

For a 7-like LSP, i.e. a x\ dominated by its photino component, the q to q + 7 coupling is 
proportional to the q electric charge and the q partial width reduces to 

1 / M-\ 2 

iW^-A^i-^) . (7) 

In such a case, more than 90% of the q — > qXi decays will involve the Xv A similar partial 
branching ratio holds for a fZ-like LSP with a relatively large Z component (e.g. in the H 
region close to the Z region in Fig. [|b). For a Z-like LSP, this branching ratio reduces to 
20% < B < 80%. Decays involving the LSP are negligible only in the H domain extending to 
negative /x's adjacent to the 7 domain (Fig. ||b). 
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(Almost) only the q~i are allowed by gauge symmetry to decay into q'xt- This is because 
the SU(2)l symmetry which implies in the SM that the right handed fermions do not couple 
to the W boson also forbids a coupling of q R to the W. The q R decays involving the chargino is 
only possible through the H + component of the x + i n which case the coupling is proportional 
to the q' mass. Hence the decay q R — > q'xt is strongly suppressed for a q R of the first or second 
generation. The partial width of the q — > qxf decay is obtained from (Eh with the interchange 



M..+ and with 



.4 



gVi 



V2 



B 



-gM q ,u l2 



(8) 



2M W cos (3 

The regions of the M 2 vs /x plane where the u decays involving a chargino dominate are shown 
in Fig. f|. In most of the parameter space, the Ul squarks will mainly undergo a decay involving 



Figure 4: Dominant gauge decay of a 
150 GeV ul squark. 
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a chargino if kinematically allowed. The mass dependence of the xt states on the \x parameter 
is shown in Fig. |^a for fixed M 2 and tan f3. 



2.4 Decays of the LSP 



In $ p SUSY models with \' x - k ^ 0, the LSP will undergo one of the following decays : Xi vd k d, 
Xi —> e + Ujdk or x\ ~^ e'Ujd^. Representative diagrams of such decays are given in Fig. || The 
relevant matrix elements for these decays can be found in fL2| . They depend on the coupling A', 
but also on the parameters M 2 , \x and tan/?. This dependence is illustrated in Fig. |^a for the 
LSP decay Xi ~^ e ± qq'. Such decay modes are seen to be dominant (63% < Br < 88%) if the Xi 
is 7-like in which case both the "right" and the "wrong" sign lepton (compared to the incident 
beam) are equally probable. This leads to largely background free striking signatures for lepton 
number violation. The latter will dominate if the Xi is Z-\ike. A H-like Xi w iU most probably 
be long lived and escape detection since its coupling to fermion-sfermion pairs is proportional 
to the fermion mass fL5f. This is illustrated in Fig. |6]b, which shows the flight distance ctq of 
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Figure 5: Example diagrams of the LSP decays Xi ~^ Iqq' involving afy p Yukawa coupling. 

the x? i n the plane (M 2 ,yu) for A' = 0.1. The cr exceeds 1 m in most of the if-like domain 
surrounding the singularity at /i = where M x o = at tree level. Hence processes involving a 
higgsino-like Xi will De affected by an imbalance in transverse momenta. 




Figure 6: (a) LSP (x\) decay branching ratio into charged leptons (i.e. e* 1 + jets), as a 
function of \i and M 2 for sfermion masses Mj = 150 GeV and tan/5 = 1; (b) log cr (m) of 
the LSP with X' = 0.1, the LSP mass is vanishingly small around \x = and along the ridge at 
large \i + M 2 . 



2.5 Decays of Charginos 



R-parity conserved x + decays into a x° an d two matter fermions, have been investigated in 



detail in |L6|, where the relevant matrix elements can be found. New decay modes of the x H 



into e + + dj + d k or v e + Uj + d k are allowed by the $ p couplings A' lj7c as illustrated in Fig. [7|. 
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Figure 7: Representative diagrams for the xf decays xt — > Iqq' involving aR\ p Yukawa coupling. 

The branching ratio of the xf m t° these R\ p modes is obtained using the partial widths 
calculated from the relevant matrix elements. Considering for instance the process x + ~^ e+ + 
dj + dk (Fig. 0a) in the following notation xt(k) ~> e+ (0 + ^(^2), and using Mandelstam 

variables s = (qi + q 2 ) 2 = (k — I) 2 , t = (k — gi) 2 = (I + (fc) 2 and u = (k — q 2 ) 2 = {I + gi) 2 , the 
squared matrix element can be written as : 

where the propagators R and D and the interference term / are : 

R(s) = s-ml L , (10) 
D(t) = t-m 2 L , (11) 
I(s,t,u) = s(M 2 -s)-u(M 2 -u) + t(M 2 -t). (12) 

The matrix element corresponding to the process xt v e + u + d (Fig. ^b) is deduced from 
the previous one with the following substitutions: e^v, d^-u, d^d and |Vn| 2 — > \Un\ 2 . 
The corresponding partial width is obtained by integrating over phase space as : 

r =/ / tt* ~m J M\ 2 dsdt (13) 

J s=0 Jt= M 3 + 256tt 3 ' 1 v ; 

The R\p decays of the xf wm mainly dominate over MSSM decays as soon as A' is not too small, 
as can be seen in Fig. [8[ For Xi — 7> $p decays of the chargino dominate over MSSM modes 
for coupling values above ~ 0.25, which is typically HERA's sensitivity limit with current 
luminosity. 



2.6 Classification of Final States 

Taking into account the dependence on the nature of the xf the possible decay chains of the 
ul and dn squarks can be classified into eight distinguishable event topologies listed in tables |] 
and and labelled SI to S8. The SI and S2 classes cover $ p squark decays. The S3 and S4 
classes are squark gauge decay topologies not accompanied by escaping transverse momenta 
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Figure 8: R\ p chargino decay branch- 
ing ratio as a function of A' m and fi, for 
M 2 = 80 GeV, tan/? = 1 and sfermions 
masses = 150 GeV; the hatched domain 
corresponds to fi values for which M(xf ) < 
M{ X l). 
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Channel 



X? 
nature 



Decay processes 



Signature 



SI 



Q 



High P ± e + + 1 jet 



S2 



H 



d R 
Q 



q 



x? 



Missing P± + 1 jet 



S3 



1,Z 



q 



A' 

d xt 



e + q'q" 



e + dd 

A' 



High Pj_ e + 
+ multiple jets 



e + q'q" 



q q 



S4 



1,Z 



1,Z 



Q Xi 

A' 



d Xi 



e q'q" 



A' 



e q'q" 



High P± e~ 

(i.e. wrong sign lepton) 
+ multiple jets 



q q 



Table 2: Squark decays in Ijt v SUSY classified per distinguishable event topologies (PART I). 
The dominant component of the Xi f or which a given decay chain is relevant is given in the 
second column. The list of processes contributing to a given event topology is here representative 
but not exhaustive. 



while those with large _Fj_ are covered by classes S5 to S8. 
has been developed and discussed in detail in §, pi] . 



A set of event selection cuts 
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Channel 



X? 
nature 



Decay processes 



Signature 



S5 



1,Z 
1,Z 

H 



Ul 



q x°i 

A' 

d Xi 



d Xi 

A' 

d Xi 



vq'q' 



A' 

U q q' 
vud 

^ q q' 



vq'q' 



Missing P± 
+ multiple jets 



d xi 



H 



W + X \ 
^1+ v 



High P± e + or fi + 
+ missing P± + 1 jet 



S7 



1,Z 



d xi 



W + X \ 

A' 



e ± q'q" 



1+ v 



High P ± e ± 
+ high P± e + or fi^ 
+ missing P± 
+ multiple jets 



S8 



1,Z 



u L 



d xi 



o 



W + x\ 
x 



uq'q' 



1+ v 



High P± e + or /x" 1 
+ missing P± 
+ multiple jets 



Table 3: Squark decays in$ p SUSY classified per distinguishable event topologies (Part II). As 
in table^, the list of processes given here is not exhaustive, e.g. the gauge decays xi 
and xi —* XiQQ' ma V a ^ so proceed via a virtual sfermion. 



For SI and S3 (or S4), the DIS NC background is strongly suppressed by requiring a high 
Pj_ e ± found at high y e , where y e = 1/2(1 + cos 9*) and 6* e is the electron angle in the e — q CM 
frame. The uniform decay of the scalar particle in the CM frame leads to a flat y e spectrum for 
SI and one shifted towards largest y e for S3. This is in contrast to the l/y^ spectrum expected 
for the DIS NC background at fixed quark momentum fraction x. For S3 the HI analysis [[l 1 



has been improved JT7], using #*'s computed for the scattered electron and for the highest P± 
jet found in the azimuthal hemisphere opposite to the electron, and cutting on V — He + yj e t 
as shown in Fig. [| Good signal detection efficiencies are obtained in each of these classes, 
reaching ~ 70% for SI and up to ~ 60% depending on M x o for S3. 

The S4 topology with a wrong sign lepton in the final state is quasi-background free. Event 
candidates in classes S2 and S5 to S8 have a large ft±. Classes S2 and S5 suffer from DIS CC 
background and from tails of photoproduction background. The S6 to S8 topologies have one 
or many leptons in the final states and are thus quasi-background free. Typical signal detection 
efficiencies ]llj reach ~ 30% — > 80% in these channels. 
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Figure 9: Distribution of the variable 
for neutral current DIS processes, and for 
a simulation of 75 GeV squarks undergoing 
gauge decays involving 20 GeV neutralinos; 
the vertical line is the cut used in the HI 
analysis / ff/|/ . 
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The relative contributions of the squark $ p and gauge decays are shown in Fig. |10|a. Gauge 
decays are seen to dominate through most of the accessible mass range. Only large Yukawa 
couplings can be probed at largest masses and thence $ p decays dominate. The shape of the 
curves in Fig. pTD| a is only distorted at lowish mass (e.g. Mq < 75 GeV) when convoluting with 
signal detection efficiencies. The measurement of the relative branching ratio in 53 and S4 
in case of a discovery, could be used to constrain the Xi LSP nature in the MSSM parameter 
space as seen in Fig. [I0|b. 

It is interesting to note that the final state classification discussed here should not be 
dramatically affected when relaxing the hypothesis of section [2.1| , e.g. in models where the g 
are lighter than the q, or where the LSP is the xt- 




Figure 10: (a) Squark decay branching ratio as a function of squark mass integrated over three 
distinct set of event topologies for tan/3 = 1 and My = 40 GeV. (b) Ratio Bsn/Bsz of the 
squark "gauge" decay branching ratios without ^± involving the like (S3) and unlike (S4) sign 
lepton viewed in the M 2 versus fi plane; the plot is obtained for Mq = 150 GeV at the expected 
limit of X' coupling sensitivity for an integrated HERA luminosity of 100 pb -1 . 
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Assuming M g < Mq, the decay q — > q + g will generally dominate. If the g is the LSP, the q 
decay will be followed by the $ p decay g — > q + q' + e ± or g — » q + q + v. In such a case, possible 
final states contain several jets and either one electron or These topologies correspond 
to channels S3 and S5, previously considered. If Mq > Mg, with the LSP being the lightest 
neutralino, the g arising from squark decay will undergo g — > q + q, the latter squark being 
off-shell. Possible final states are similar to those listed above, but more jets would be expected. 
Assuming now that the LSP is the xt ( see the relevant MSSM parameters in Fig. |3|b), a new 
event topology would only emerge for a relatively stable xt which could behave as a "heavy 
muon". However, the time of flight of the xt> obtained from the integration (|T^) over phase 
space, reads as : 

4vr 1 1 . 2x / M s \ 4 1 

(8x64vr 2 )U^ — - (14) 



g 2 |Vii| 2 A /2 v ' \ M v+ j M v+ 



which numerically leads to : 



From this formula one obtains that the relevant parameter space for the xt to decay outside the 
detector (> lm), is already excluded from the intrinsic Z° width measurement at CERN |i 



3 Results for the Mass- Coupling Reach of HERA 

In the absence of a significant deviation from the SM expectations, exclusion limits for the 
Yukawa couplings A' ljfc as a function of mass can be derived, showing the domain HERA could 
probe in the near future. Results are shown for X' 1 j 1 in Fig. |Tl| at 95% confidence level (CL), for 
integrated luminosities L = lOOpb -1 and L = 500pb _1 . These have been obtained assuming 
a 40 GeV 7-like x?> an d combining all contributing channels. For C = 500pb _1 , the existence 
of first generation squarks with fi p Yukawa coupling X'^ could be excluded for masses up to 
~ 270 GeV for coupling strengths Af 1;L /47r > a em . 

From the analysis of the X' 1 j 1 case involving the cIr and ul squarks, limits can be deduced 
on the A' ljfc by folding in the proper parton densities. Such limits are given in Table § at 
Mq = 150 GeV and for an integrated luminosity of 500pb~\ Also quoted in this table are 
the most severe existing indirect limits for each of these couplings. The most stringent concern 
couplings X[j k with j — k and come either from the non-observation of neutrinoless double-beta 
decay (j = k = 1) or from constaints on the v e mass (j = k = 2, 3). The limit from (3f3Qv decay 
depends on the gluino mass and is given here for M g = 500 GeV. 

By the time HERA reaches high luminosity running conditions, new direct limits (or a 
discovery !) from other colliders will have further constrained the possible squark masses 
and SUSY parameters. In e + e~ collisions, the direct squark pair production process does not 
violate R-parity and LEP2 should directly probe squark masses up to \fs/2, i.e. ~ 90 GeV. 
In pp collisions, squarks can be produced in pair or in association with gluinos. No complete 
analysis in the ^? P -SUSY framework has been performed yet with existing TEVATRON data. 
Nevertheless, q decay topologies similar to those described here have been explored by DO 



and CDF [26] in scalar leptoquark or MSSM searches. From these and from di-lepton data [27] , 
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Figure 11: (a) Exclusion upper limits at 95% CL on the X' l j l coupling as a function of squark 
mass which could be reached with e + p collisions at HERA (^fs ~ 300 GeV) for integrated 
luminosities of L = 100 pb _1 (dark shaded area) and 500 pb _1 (shaded); (b) Regions of the 
M-2 versus /j, plane excluded for L = 100 pb _1 and for couplings X' 1 j 1 equal or smaller to the 
exclusion upper limit at Mg = 150 GeV. 

Table 4: Exclusion upper limits at 95% CL on the coupling A' ljfc for Mg = 150 GeV and M x o = 
40 GeV together with best existing indirect limits. The indirect limits have been scaled from 
those found in the cited references to Mg = 150 GeV and 95% CL. 



HERA sensitivity Indirect limits 





7-like Xi 


Z-like x \ 


Value [Ref. 


Nature of the process 




0.008 


0.023 


0.003 


y 


(3(30v decay 


•^112 


0.020 


0.057 


0.05 


m 


CC-universality 


^113 


0.026 


0.072 


0.05 


m 


CC-universality 


^121 


0.008 


0.023 


0.09 


21 


Atomic Parity Viol. 


\/ 

A 122 


0.027 


0.077 


0.04 


22j 


i/ e -mass 


\' 

A 123 


0.043 


0.012 


0.5 


m 


D + — > K decays 


^131 


0.007 


0.024 


0.09 


21 


Atomic Parity Viol. 


\' 

A 132 


0.027 


0.091 


0.77 


24 


RfP 


X' 

A 133 


0.068 


0.230 


0.0015 


22| 


z/ e -mass 



one can infer that the range 200 — > 300 GeV of ^? P -SUSY squark masses will most probably 
be not fully excluded by TEVATRON data for an integrated luminosity of ~ 100 pb _1 , thus 
leaving open a discovery window at HERA in the hypothesis Mg 3> Mg. 

If the presence of two simultaneously non-vanishing Yukawa couplings (e.g. A' ljfc and A- jfc 
with j ^ 1), resonant q production at HERA can be directly followed by a lepton flavor violation 
(LFV) decay leading to /x + jet or r + jet signatures. Relevant analysis with existing data have 
been performed by the HI |13[ and ZEUS |28| collaborations and limits comparable to the best 
existing indirect LFV limits have been derived in the context of $ P -SUSY for a pure 7 like LSP. 
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A new range of possible coupling products could be probed with increasing luminosity 



29 



4 Summary and Conclusions 

The HERA potential for i?-parity violating supersymmetry searches was studied. Direct res- 
onant production of squarks through nine new Yukawa couplings A^- fe is possible up to the 
kinematical limit of ~ 300 GeV. 

Supersymmetric partners qi of left handed w-like squarks are produced preferentially in e + p 
collisions and most favourably via A' lj:L . In contrast, e~p collisions mainly produce partners qn 
of right-handed d-\ike quarks and most favourably via \' llk . S quark decays via a A' coupling into 
/ + q final states dominate only at largest accessible masses, while elsewhere squarks undergo 
mainly gauge decays into a quark and a gaugino-higgsino. The qu decays involve a neutralino 
X° while q~i decays dominantly proceed via a chargino x + i n a large portion of the MSSM 
parameter space. The x' s > including the LSP, are generally unstable and their decay chain 
involves the A' ljfc coupling. 

In total, eight classes of event topologies are identified for i?-parity and gauge decays of 
squarks, with single or multi-leptons final states accompanied or not by missing transverse 
momenta. A good experimental sensitivity is expected in each of these classes. Thus, for an 
integrated luminosity of 500 pb _1 , squarks can be searched for Yukawa couplings smaller than 
the electromagnetic coupling up to masses of < 270 GeV, almost independently of the specific 
choice of MSSM parameter values. Coupling values below the most stringent indirect limits 
can be probed at Mg = 150 GeV for seven out of the nine possible A' ljA couplings. 
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Appendix: Gaugino-higgsino mixing 



Detailed expressions (and conventions) used for the mass mixing in the gaugino-higgsino 
sector are presented here for completeness. 

Neutralino mass mixing: 

Mass terms of the Lagrangian describing SU(2)l X U(1)y neutral gauginos and higgsinos can 
be written as : 

C m = -l/2(^ ) T y^° + h.c. 
where the neutralino mass matrix in the basis ip® = (— iA, —iW 3 , H®, ifg) is given by 

/ 



Y 



Mi 
M 2 

—mz sin 9w cos (3 mz cos 9w cos (3 

V mz sin 9w sin (3 —mz cos 9w sin (3 



-mz sin 9w cos (3 m z sin 9 W sin (3 \ 

sin 9w cos /3 -mz cos 9w sin /3 


-H / 



(16) 



(17) 



The number of free MSSM parameters is reduced by using a GUT inspired relation between 
soft-breaking terms Mi and M 2 , Mi = | tan 2 6> W M 2 . 

Neutralinos correspond to the mass eigenstates and are defined as Xi = N^ipj, with iVy 
being the unitary matrix which diagonalize Y . Finally, we make use of the matrix N', which 
diagonalizes the neutralino mass matrix expressed in the basis (7, Z) instead of (A, W3) : 
iVji = Nj! cos 9 W + N j2 sin 9 W , N' j2 = -N n sin 9 W + N j2 cos 9 W , N' 



"j3 



N j3 and N' j4 



N. 



J'4- 



Chargino mass mixing: 

The Lagrangian mass terms for winos and charged higgsinos are written as : 



*— "r 



X T 
X 



+ h.c 



(18) 



where : 



X 



M 2 m w \/2sm(3 
mvi/v^2cos/5 /x 



(19) 



and ipj~ = (— iW + , H^), ipj = (—iW ,H\ )- This mass matrix is diagonalized using two (2, 2) 



unitary matrices U and V |L6| : = V lJ ipj and Xi = U lJ ipj . Masses for these eigenstates 



are easily derived from the above X matrix : 



w 



M 2 2 = |(M| + /i 2 + 2m : 

=F [(Mf - fi 2 ) 2 + Am A w cos 2 2(3 
+4m 2 w (M 2 + /i 2 + 2M 2y usin2/3)] 1/2 



(20) 
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